T fd E HH T 


ChinaxX iv HA 


FE 


PXA 


F 


ARID ZONE RESEARCH 


doi; 10. 13866/j. azr. 2019.02.22 


基于 3S 技术 的 天 山 历 史 云 杉林 空间 分 布 的 提取 


邢 菲 ， 李 虎 ， 李 建 贵 ， 


aK, AER, RAW 


(1. 新疆 农业 大 学 草 业 与 环境 科学 学 院 ,新 疆 乌鲁木齐 ”830052; 2. 潍 州 学 院 地 理 信息 与 旅游 学 院 ,安徽 滁州 239000 ; 
3. 新 疆 农 业 大 学 林业 研究 所 ,新 疆 乌鲁木齐 ”830062; 4. 新 疆 师范 大 学 地 理科 学 与 旅游 学 院 ,新 疆 乌鲁木齐 ”830054 ) 


摘 要 : 运用 遥感 手段 结合 历史 时 期 遥感 影像 数据 ,以 天 山 云 杉 ( Picea schrenkiana var. tianschanica ) 林 生境 特征 为 


的 同时 提高 运行 速度 。 从 整体 来 看 , 利 


固定 因子 ,结合 植被 指数 分 析 ` 地 形 因子 分 析 Eo ot 
林 的 空间 分 布 信息 ,从 而 为 历史 资料 缺失 情境 下 的 天 然 林 资源 保护 工程 实施 效益 评价 提供 支持 。 研 究 表明 :GD 将 
天 山 云 杉林 的 年 龄 特征 作为 固定 因子 ,以 现状 年 高 空间 分 辩 率 遥感 影像 及 森林 资源 二 类 调查 数据 作为 本 底 资 料 ， 
在 面向 对 象 分 类 方法 支持 下 可 以 很 好 的 从 历史 时 期 的 遥感 影像 中 提取 出 天 山 云 杉林 的 历史 空间 分 布 信息 ,提取 精 
度 可 达 93.3% 52) 在 植被 指数 因子 中 ,NDVI 对 天 山 云 杉林 指示 性 最 好 ,并 确定 用 于 天 山 云 杉林 提取 的 最 佳 NDVI 
值 域 为 [0.35 ,0.8] ;@) 地 形 因子 及 主 成 分 分 析 方 法 可 以 大 大 压 双 
用 历史 时 期 遥感 影像 并 结合 天 山 云 杉林 的 生境 特征 ,可 以 很 好 的 提取 出 历 


向 对 象 的 决策 树 分 类 的 方法 提取 历史 时 期 天 山 云 杉 


影像 的 元 余 信 息 , 在 提升 云 杉 林 信 息 提取 的 精度 


史 时 期 云 杉林 空间 分 布 信息 ,从 而 为 资料 缺失 情境 下 的 森林 资源 管理 措施 制定 及 应 对 气候 变化 提供 数据 支持 。 
关键 词 : 历史 遥感 影像 ; 决策 树 分 类 ; 天 山 云 杉林 ; 空间 分 布 ; 量 康 林场 


森林 既是 历史 自然 环境 的 反映 者 ,也 是 生物 圈 
中 最 富有 生机 和 巨大 生产 力 的 生态 系统 ""。 位 于 
干旱 区 的 森林 更 是 范 漠 中 的 “ 湿 岛 ”", 对 于 改善 周边 
局 地 气候 ,保持 土壤 和 涵养 水 源 具 有 显著 的 作用 。 
天 山 云 杉 是 亚洲 中 部 山地 的 特有 树种 ,在 我 国 仅见 


查 资料 ,将 云 杉林 分 布 信息 与 空间 分 析 技 术 相 结 
BOM ;并 综合 考虑 云 杉林 的 生理 特征 地形、 地 
貌 特征 ;利用 植被 在 遥感 影像 上 的 光谱 信息 特 
PEUS ;以 历史 遥感 影像 为 数据 源 ,在 RS 软件 平台 支 
持 下 利用 面向 对 象 分 类 模型 ,提取 天 山 云 杉林 历史 


于 新 疆 , 并 且 是 新 疆 山 地 森林 中 分 布 最 广 .蓄积 量 最 
大 的 树种 中 。 天 山 云 杉林 是 天 山 山地 最 主要 的 地 
带 性 植被 中 ,以 天 山 云 杉 为 主体 的 天 山 山地 森林 生 
态 系统 对 水 源 的 涵养 起 着 至 关 重 要 的 作用 '"。 天 
山 云 杉 在 新 疆 主 要 分 布 于 天 山北 贰 的 西 坡 、 北 坡 和 
西北 坡 ;天 山南 苑 分 布 较 少 ;在 昆仑 山西 部 北 坡 也 有 
少量 分 布 。 

目前 ,关于 天 山 “- 及 天 山 云 杉 的 种 群 结构 特 
征 "-" 生物 量 及 变化 规律 ~- 碳 储 量规 
律 "和 生境 特征 “等 已 有 较 多 研究 。 但 天 山 
云 杉 林 的 历史 分 布 情况 如 何 , 在 相关 历史 资料 遗失 
的 情况 下 云 杉林 的 分 布 如 何 进行 调查 ?基于 此 ,本 
文 在 缺乏 历史 调查 数据 的 背景 下 ,以 天 山北 坡 中 段 
息 康 林 区 天 山 云 杉林 为 研究 对 象 ;结合 现 有 林业 调 
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空间 分 布 信息 ,在 为 利用 遥感 技术 提取 森林 信息 提 
供 了 一 种 新 的 思路 的 同时 ,为 森林 资源 可 持续 经 营 
管理 和 应 对 全 球 气候 变化 提供 数据 支持 。 


1 研究 区 概况 


WEE XK VE FE LL P BAL FE 88°00" ~ 88°42'E， 
43°38 ~44°05'N 的 阜 康 林场 , 量 康 林场 始 建 于 1956 
年 3 月 ,总 面积 92 750. 86 hm (图 1) 。 原 称 阜 康 县 
四 工 森林 经 营 所 ,林场 经 历 6 次 变革 后 至 2013 年 转 
为 全 额 事业 单位 ,归属 和 旦 康 市 管辖 。 多 次 的 隶属 关 
系 变更 导致 林场 内 历史 数据 资料 丢失 严重 ,无 法 满 
足 现 有 森林 经 营 规划 及 林业 保护 工程 工作 的 需要 。 
研究 区 东 与 吉木 萨 尔 林场 相连 , 西 与 米泉 林场 相 邻 ， 
南 为 东 天 山 主 峰 博 格 达 峰 , 北 接 准噶尔 盆地 。 研究 
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图 1 研究 区 示意 图 
Fig.1 Sketch of the study area 


DX PY BRAKE BEE FR LL zs 42 9 CY LL SE a PEEP 
叶 林 ,森林 带 出 现在 海拔 1 500 ~2 700 m 的 坡 面 上 ， 
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面积 达 15 841.45 hm (2016 年 数据 ) 。 研 究 区 属于 
温带 大 陆 性 干旱 、 半 干旱 气候 ,年 平均 气温 1.2 C, 
年 降水 量 400 mm ,降水 年 内 分 配 不 均 ,主要 分 布 在 
夏季 和 冬季 。 研 究 区 土壤 分 布 呈 明显 的 垂直 地 带 
性 ,由 低 到 高 依次 为 :山地 栗 钙 土山 地 灰 褐色 和 森林 
土 .高山 草 多 土 和 高 山石 质 带 。 


2 数据 与 方法 


2.1 数据 源 与 预 处 理 
本 研究 所 用 的 数据 有 :遥感 影像 数据 .森林 资源 
二 类 调查 数据 研究 区 数字 高 程 模型 (DEM ) 、 植 被 
首 数 坡度 和 坡 向 数据 。 数 据 的 来 源 、 参 数 及 处 理 方 
法 见 表 1。 


表 1 数据 源 、 参 数 及 处 理 方 法 


Tab.1 Data sources, parameters and processing methods 


数据 类 型 SRR 数据 来 源 及 处 理 
rer O. 8m 高 分 2 号 过 感 影像 数据 (2016 年 ) 来 源 于 中 国资 源 卫星 应 用 中 心 网 站 , 利 
erie ene 0-8 m2016 年 。 GIS 和 RS 平台 进行 影像 的 校正 ,拼接 及 研究 区 范围 裁 前 

E 30 my1904 ip 30m SPER Landsat TM 数据 来 源 于 USGS 网 站 ,利用 GIS 和 RS 平台 进行 影像 的 校 


正 、 拼 接 及 研究 区 范围 裁剪 


2016 年 森林 资源 二 类 调查 数据 10 000/2016 年 ” 阜 康 林场 提供 ,利用 GIS 平 台 对 数据 进行 规范 化 处 理 ,并 提取 云 杉 林 信 息 
数字 高 程 模型 (DEM ) 30 m 来 源 于 地 理 空间 数据 云 网 站 ,利用 GIS 平台 进行 拼接 和 研究 区 范围 裁剪 
坡度 30 m 利用 GIS 平台 根据 研究 区 DEM 数据 提取 
坡 向 30 m 利用 GIS 平台 根据 研究 区 DEM 数据 提取 
植被 指数 ( NDVI) 30 m/1994 年 。 利用 RS 软件 平台 根据 1994 年 遥感 影像 提取 


2.2 采样 点 选择 

以 研究 区 2016 年 国产 高 分 2 号 0.8 m 分 辨 率 
遥感 影像 和 从 2016 年 开展 的 森林 资源 二 类 调查 数 
据 中 提取 的 云 杉林 为 本 底 资 料 , 根 据 云 杉林 生长 慢 
的 特点 ,剔除 龄 组 为 幼 龄 林 的 云 杉林 ,保留 林 龄 为 中 
龄 林 、 近 熟 林 、 成 熟 林 和 过 熟 林 的 云 杉林 以 保证 云 杉 
林 在 近 30 a 内 不 发 生变 化 ,并 在 每 个 林 班 内 随机 选 
择 1~2 个 小 班 的 中 心 点 作为 样 点 数据 (不 区 分 龄 
级 ) , 共 选 择 100 个 样 点 (图 1) , 样 点 中 主要 包括 经 、 
纬度 信息 ,树龄 信息 ,树种 信息 。 其 中 经 纬度 信息 主 
要 用 于 与 地 形 数据 及 植被 指数 到 加 时 位 置 匹 配 , 树 
龄 和 树种 信息 作为 云 杉林 生长 的 固定 因子 ,以 保证 
该 处 在 历史 上 也 是 云 杉林 不 变 。 并 依据 分 类 样 点 与 
仿 验 样 点 7:3 的 原则 ,选择 70 个 样 点 作为 分 类 样 点 
用 于 天 山 云 杉林 信息 提取 ,30 个 样 点 作为 检验 样 点 
进行 精度 验证 。 
2.3 植被 指数 提取 

遥感 图 像 上 的 植被 信息 主要 通过 绿色 植物 叶子 


和 植被 冠 层 的 光谱 特性 及 其 差异 .变化 反映 。 从 
遥感 数据 中 获取 的 植被 指数 已 经 被 广泛 用 来 定性 和 
定量 评价 植被 覆盖 度 及 其 生长 活力 ”-”。 本 文选 
取 3 种 常用 的 植被 指数 :比值 植被 指数 (RW) 、 归 一 
化 植被 指数 (NDV1) 和 增强 型 植被 指数 (EVI) 作 为 
天 山 云 杉林 提取 的 候选 因子 ,通过 对 比 天 山 云 杉林 
对 不 同 植被 指数 的 响应 程度 ,最 终 选 取 一 种 植被 指 
数 作为 云 杉林 提取 的 参数 因子 。3 种 植被 指数 的 提 
取 方 法 如 下 : 
RVI =P e Prix (1) 
NDVI = (Prit ~ Prea )/ (Puit + prea) (2) 
EVI=G[ (Prr = Prea )/ (Pair + APrea + Asps +L) (3) 
式 中 : pi 为 红 波 段 的 地 表 反射 率 ;pi 为 近 红 外 波段 
FAIS AS ST AR spy, Fy WE UG Be HY) St Ae I Bs GA, A, 
A L Ay HEX 2.5 6.0.7.5 和 1.0。 
2.4 地 形 因子 提取 
参照 杨 听 、 刘 娇 等 “ ”的 理论 与 方法 ,在 GIS 
软件 平台 支持 下 提取 研究 区 内 坡 向 和 坡度 信息 。 并 
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WS ABS AEF 3S 技术 的 天 山 历 史 云 杉林 空间 分 布 的 提取 


(al) DEM 


(a2) AHESDEMB In 


图 2 研究 区 地 形 因子 


Fig.2 Topographic factors in the study area 


以 北方 向 为 起 点 , 顺 时 针 每 隔 4S RFE Tel oo AE A 
东 、 东 、 东 南 、 南 、 南 西 . 西 .西北 8 个 方向 ;根据 坡度 
值 <5°、5° ~ 14°、15° ~ 24°、25° ~ 34°、35° ~ 44° 
和 >45° 将 坡度 分 为 平 \ 缓 .和 斜 \ 陡 、 急 、 险 6 个 等 级 ， 
研究 区 DEM \ 坡 度 \ 坡 向 如 图 2al ,2b1 ,2cl。 
2.5 主 成 分 变换 

主 成 分 分 析 ( 又 称 K-L 变换 ) 是 一 种 除去 波段 
之 间 的 多 余 信息 ,将 多 波段 的 图 像 信息 压缩 到 比 原 
波段 更 有 效 的 少数 几 个 转换 波段 的 方法 E RS 
软件 平台 的 支持 下 , 主 成 分 变换 操作 简单 .处理 速度 
快 ,在 遥感 数据 信息 提取 中 已 经 取得 了 较 好 的 应 用 
BORE 。 通 过 主 成 分 变换 ,可 以 去 除 噪声 对 地 物 
信息 的 影响 ,将 影像 的 主要 信息 压缩 到 第 一 主 成 分 
中 ,增强 地 物 的 细部 特征 ,更 有 利于 对 地 物 的 分 析 与 
提取 。 本 研究 在 RS 软件 平台 的 支持 下 ,对 人 研究 区 
影像 进行 主 成 分 分 析 ,以 降低 影像 噪声 的 影响 ,从 而 


突出 天 山 云 杉林 信息 。 
3 结果 与 分 析 


3.1 研究 区 植被 指数 

根据 植被 指数 提取 方法 ,在 RS 软件 支持 下 根 
据 历史 遥感 影像 (1994 年 ) 提取 得 到 研究 区 EVI, 
NDVI 和 RVI 3 种 植被 指数 (图 3 ) ,并 将 70 个 采样 点 
分 别 与 3 种 植被 指数 进行 琶 加 分 析 , 统 计 各 个 样 点 
所 在 位 置 植被 指数 的 最 大 值 . 最 小 值 平 均值 标准 
差 , 并 分 别 计算 各 采样 点 与 EVINDVI 和 RVI 的 相 
关系 数 , 结 果 如 表 2。 

从 结果 的 最 大 值 最 小 值 来 看 , EVI NDVI 和 RVI 
3 种 植被 指数 均 在 0 ~1 范围 内 ,从 平均 值 来 看 ， 
NDVI 的 均值 为 0.75 ,而 EVI 和 RVI 的 均值 分 别 为 
0.21 和 0. 15; 从 云 杉 林 采 样 点 与 3 种 不 同 植被 指数 
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TE EVI 为 增强 植被 指数 ,NDVI 为 归 一 化 植被 指数 ,RV 为 比值 植被 指数 。 


AL 3 植被 指数 提取 结果 (1994 年 ) 


Fig.3 Extracted results of vegetation indexes in 1994 


表 2 云 杉 林 与 植被 指数 又 加 分 析 结 果 
Tab.2 Coupled results of vegetation indexes 
植被 指数 类 型 ”最 小 值 ” 最 大 值 ” 平 均值 ”标准 差 相关 系数 
EVI 0. 14 0. 44 0. 21 0. 06 0. 36 
NDVI 0. 03 0. 87 0.75 0. 05 0. 45 
RVI 0. 07 0. 23 0.15 0.04 -0.02 
TE EVI 为 增强 植被 指数 ,MD 人 为 归 一 化 植被 指数 ,RVI 为 比值 植被 


间 的 相关 关系 来 看 , NDVI > EVI > RVI, HRERS 
NDVI 相关 系数 (0. 45) 要 明显 大 于 采样 点 与 EVI 
(0. 36) A RVI( -0. 02) 。 综 合 以 上 分 析 结 果 ,并 将 
NDVI 用 于 提取 实验 ,最 终 得 出 NDVI 在 值 域 [0. 35, 
0.8] 时 天 山 云 杉林 的 提取 效果 最 好 。 
3.2 研究 区 地 形 因子 

根据 研究 区 海拔 坡度 和 坡 向 提取 结果 ,利用 
GIS 平台 的 空间 分 析 功 能 将 研究 区 森林 资源 二 类 调 
查 数据 与 地 形 因子 (图 2a2 ,2b2 ,2c2 ) ap HET IM 
分 析 得 到 结果 如 图 4。 

从 研究 区 天 山 云 杉林 分 布 海拔 来 看 , 有 
96. 24% 的 天 山 云 杉林 分 布 在 海拔 1 600 ~2 800 m 
范围 内 , 且 在 海拔 1 000 ~2 200 m 范围 内 天 山 云 杉 


林 的 分 布 随 海 拔高 度 的 上 升 而 增加 ,并 在 海拔 2 200 
m 达到 峰值 ,而 后 随 着 海拔 高 度 的 上 升 而 减少 。 从 
云 杉林 分 布 坡度 来 看 ,在 陡坡 分 布 范围 最 广 , 占 总 分 
布 范围 的 36. 91% , 其 次 是 斜坡 和 和 急 坡 ,分别 占 
26.34% 和 21. 03% ,缓坡 分 布 较 少 , 占 11. 36% , 共 
Th 95. 65% 的 天 山 云 杉林 分 布 于 缓坡 BEDE BER A 
急 坡 。 从 天 山 云 杉林 坡 向 分 布 来 看 ,共计 85. 49% 
的 天 山 云 杉林 分 布 于 东 坡 、 东 北 坡 、 西 北 坡 和 北 坡 ， 
其 中 北 坡 和 东北 坡 分 布 最 广 ,分 别 占 26. 72% 和 
24.86% .结合 研 究 区 历史 遥感 影像 及 地 形 因子 分 析 
结果 ,确定 海拔 区 间 [1 550,2 890 ] ,坡度 区 间 [5 
50] , 坡 向 区 间 [0 ,135 ] 和 [225 ,360 ] 作 为 天 山 云 杉 
林 历 史 空 间 分 布 提取 因子 。 
3.3” 主 成 分 分 析 
根据 主 成 分 计算 方法 ,在 RS 平台 支持 下 计算 
影像 各 个 主 成 分 ( 表 3 ) ,其 中 前 三 个 主 成 分 包含 了 
99.76% 的 信息 量 。 将 信息 量 小 而 噪声 大 的 主 成 分 
波段 舍弃 ,保留 信息 量 最 丰富 的 第 一 、 第 二 和 第 三 主 
成 分 ,并 将 采样 点 与 主 成 分 波段 进行 耦合 分 析 ( 表 
3)。 其 中 ,第 一 主 成 分 在 值 域 [ -6,40] , 占 总 体 样 
点 的 90.75% , 在 大 于 40 以 后 虽然 占 总 体 样 点 的 
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图 4 地 形 因子 分 析 


Fig.4 


Analysis on topographical factors 
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RI ， 样 点 与 第 一 第 二 及 第 三 主 成 分 耦合 结果 
Tab.3 The results of sampling sites coupled 
with PCA1 ,PCA2 and PCA3 


PCAL PCA2 PCA3 
值 域 占 比 /% 值 域 ” 占 比 /% 值 域 ” 占 比 /% 
-6~9 43.40 -1~10 56.60 -23~17 45.28 
9 ~24 33.96 10~21 30.19 -17~10 32.08 
24 ~39 13.21 21 ~32 ye) -10~3 11.32 
39 ~54 5. 66 32 ~54 3.77 -3~3 9.43 
54 ~59 ST 54 ~65 1.89 3 ~10 1.89 


9.43% ,但 分 布 离散 ,在 39 ~ 59 范围 内 只 有 5 个 样 
点 分 布 ;第 二 主 成 分 中 值 域 [ -1,32] , 占 总 体 样 点 
的 94.34% ,在 32 以 上 只 占 总 体 样 点 的 3. 66% 且 分 
布 离散 ,在 32 ~ 65 范围 内 只 有 3 个 样 点 分 布 ;第 三 
主 成 分 中 值 域 [ -23 ,3] , 占 总 体 样 点 的 98. 11% ,在 
3~ 10 范围 内 只 有 1 个 样 点 分 布 , 占 总 体 样 点 的 


} 
坡 
坡 


分 类 输入 影像 (1994 年 ) 


P 


o es 3 -6 < PCA1 <40 
a ee 0.35<NDVI<0.8 |} -1<PCA2<32 
(0° < PRIA < 135°)& 

225° < 坡 向 < 360° -23 <PCA3 <3 


{OQ 
天 山 云 杉林 提取 结果 


1.89% 。 
3.4 决策 树 分 类 与 结果 后 处 理 

决策 树 分 类 作为 一 种 监督 分 类 的 专家 分 类 方 
法 。 其 方法 灵活 直观 结构 清 晰 ,在 处 理 多 维 问题 
上 具有 传统 方法 所 不 及 的 优势 ,而 且 对 于 不 同 数据 
源 的 植被 分 类 具有 一 定 的 普 适 性 ,这 些 特点 使 得 决 
策 树 分 类 方法 在 各 种 遥感 信息 提取 中 得 到 广泛 应 
用 SFE RS 软件 平台 支持 下 ,将 植被 指数 、 海 
R .坡度 , 坡 向 及 主 成 分 影像 进行 波段 组 合作 为 分 类 
输入 影像 ,结合 研究 区 地 形 因 子 .植被 指数 和 主 成 分 
分 析 结 果 ,在 RS 软件 平台 支持 下 构建 分 类 决策 树 
规则 (图 5) 。 根 据 决策 树 分 类 规则 ,在 RS 软件 中 构 
建 决策 树 分 类 模型 并 执行 ,对 模型 提取 结果 进行 分 
类 后 处 理 得 到 天 山 云 杉 林 历 史 空 间 分 布 提取 结 
(图 5)。 


云 杉林 历史 空间 分 布 提取 结果 


云 杉林 


0 5 10km 
mm 


ALS 决策 树 分 类 规则 及 天 山 云 杉林 历史 空间 分 布 提取 结 


Fig.5 Classification rules of the decision tree and the extracted results of the history distribution of Picea schrenkiana Rupr 


表 4 决策 树 分 类 验证 


Tab.4 Verification of the decision tree classification 


耦合 结果 样 点 个 数 占 比 /% 
云 杉林 28 93.3 
非 云 杉林 2 6.7 


3.5 分 类 精度 验证 

以 30 个 检验 样 点 构建 感 兴趣 区 并 利用 监督 分 
类 中 的 最 大 似 然 法 进行 精度 评价 ,结果 显示 :最 大 似 
然 法 结果 总 体 分 类 精度 为 90. 93% , Kappa 系数 为 
0.89。 同 时 ,将 相同 的 检验 样 点 与 天 山 云 杉林 提取 
结果 进行 耦合 ,记录 每 个 检验 样 点 所 在 位 置 的 提取 
结果 并 进行 统计 ,结果 显示 :在 所 有 30 个 检验 样 点 
中 有 28 个 样 点 分 为 云 杉林 , 占 总 体检 验 样 点 的 


93.3% ;2 个 样 点 被 分 为 非 云 杉林 , 占 总 检验 样 点 的 
6.7% ;对 比 监督 分 类 结果 ,决策 树 分 类 精度 高 出 
2.37 个 百分点 ,如 表 4 所 示 。 


4 ”结论 与 讨论 


(1) 将 天 山 云 杉林 与 地 形 因子 钱 加 分 析 时 得 出 
如 下 结果 :天 山 云 杉林 分 布 海拔 范围 为 1 550 ~ 
2 890 m , 坡 向 为 5$" ~50°* 和 225° ~360° 的 东北 坡 、 西 
北 坡 和 少量 西南 坡 , 坡 度 在 5° ~ 50° 之 间 。 海 拔 范 
围 . 坡 向 略 宽 于 已 有 人 研究 成 果 ,坡度 与 其 他 学 者 的 研 
FEB AAR 。 

(2) 在 几 种 常用 的 植被 指数 中 ,对 比 EVI NDVI 
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及 RVI 与 云 杉林 的 相关 性 ,VDT7 对 天 山 云 杉林 的 信 
息 表 达 最 为 明显 , 且 用 于 天 山 云 杉林 提取 的 最 佳 
NDVI 值 域 为 [0. 35 ,0. 8] 。 

(3) 结合 云 杉林 生境 特征 ,使 用 决策 树 分 类 方 
法 对 云 杉 林 分 布 信息 提取 精度 可 达 93.3% , 比 使 用 
传统 监督 分 类 方法 精度 高 2. 37 个 百分点 ,可 以 准确 
的 提取 出 天 山 云 杉林 的 历史 空间 分 布 信息 。 

(4) 通过 分 类 精度 评价 及 其 与 传统 方法 比较 ， 
本 研究 所 使 用 方法 可 以 有 效 提 取 历 史 时 期 的 天 山 云 
杉林 分 布 信息 ,从 而 解决 历史 数据 丢失 的 问题 ,为 森 
林 经 营 管理 规划 和 应 对 气候 变化 提供 数据 支撑 。 

本 研究 使 用 的 决策 树 分 类 方法 对 云 杉林 历史 时 
期 的 空间 分 布 信息 的 提取 ,取得 了 较 好 的 结果 。 本 
研究 区 内 云 杉林 分 布 特征 与 其 他 学 者 研究 成 果 稍 有 
出 入 ,可 能 与 天 山 云 杉 自身 地 理 分 布 差异 及 使 用 的 
数据 精度 有 关 。 通 过 本 研究 使 用 的 方法 可 以 提取 出 
天 山 云 杉林 的 历史 分 布 信息 , 却 无 法 得 到 云 杉 林 的 
小 班 区 划 信 息 和 蓄积 量 信息 ,这 都 有 待 于 深入 探讨 。 
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Abstract: The spatial distribution information of Picea schrenkiana var. tianschanica in the Tianshan Mountains 
in historical period was extracted based on the vegetation index, topographic factor, principal component analysis , 
the decision tree classification method and the habitat characteristics of P. schrenkiana in the study area using the 
remote sensing methods combined with the historical remote sensing image data. So as to provide support for the 
benefit evaluation of the natural forest resources protection project under the situation of missing historical data. Re- 
sults showed that; @) the historical spatial distribution information of P. schrenkiana in the Tianshan Mountains 
could be extracted from the remote sensing images, the forest stand age of P. schrenkiana was set as a fixed factor, 
and the present remote sensing images with high spatial resolution and forest management investigation data were 
used as the background information. The accuracy of information extraction of P. schrenkiana in the study area 
could be as high as 93.3% , and the remote sensing images can be used to extract the spatial distribution informa- 
tion of P. schrenkiana in the Tianshan Mountains;@) In the vegetation index factors, the response of P. schrenki- 
ana in the Tianshan Mountains to NDVI was the most sensitive, and the best NDVI range for extracting the informa- 
tion of P. schrenkiana in the Tianshan Mountains was [0.35, 0.8]; @) Topographic factor and principal compo- 
nent analysis method could greatly compress the redundant information of image, which improved the accuracy of 
information extraction of P. schrenkiana forest and improved the running speed. On the whole, the spatial distribu- 
tion information of P. schrenkiana forest in the historical period can be well extracted by using the historical remote 
sensing images and combining with the habitat characteristics of P. schrenkiana forest in Tianshan Mountains, so as 
to provide data support for the formulation of forest resource management measures and the response to climate 
change in the context of data missing. 
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